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AN ALGORITHM FOR CONSTRUCTING CERTAIN DIFFERENTIAL 
OPERATORS IN POSITIVE CHARACTERISTIC 

ALBERTO F. BOIX*, ALESSANDRO DE STEFANI f , AND DAVIDE VANZO 


Abstract. Given a non-zero polynomial / in a polynomial ring R with coefficients in a finite field 
of prime characteristic p, we present an algorithm to compute a differential operator $ which raises 
1 // to its pth power. For some specific families of polynomials, we also study the level of such a 
differential operator <5, i.e., the least integer e such that 5 is R p -linear. In particular, we obtain a 
characterization of supersingular elliptic curves in terms of the level of the associated differential 
operator. 


1. Introduction 

Let R = k[x i, ... ,xy] be the polynomial ring over a field k, and let Vr be the ring of fc-linear 
differential operators on R. For every non-zero / £ R, the natural action of Vr on R extends 
uniquely to an action on Rf. In characteristic 0, it has been proven by Bernstein in the polynomial 
ring case (cf. [Ber72, Corollary 1.4]) that Rf has finite length as a D^-module. The minimal m such 
that Rf = Vr ■ -jr r is related to Bernstein-Sato polynomials (cf. [ILL + 07, Theorem 23.7, Definition 
23.8, and Corollary 23.9]), and there are examples in which m > 1 (e.g. [ILL + 07, Example 23.13]). 
Remarkably, in positive characteristic, not only is Rf finitely generated as a Xbj-module [Bpg95, 
Proposition 3.3], but it is generated by j (cf. [AMBL05, Theorem 3.7 and Corollary 3.8]). This is 
shown by proving the existence of a differential operator 5 £ Vr such that 5(l/f) = l// p , i.e., 
a differential operator that acts as the Frobenius homomorphism on 1//. The main result of this 
paper exhibits an algorithm that, given f £ R, produces a differential operator 5 £ Vr such that 
5(l/f) = 1 / f p . We will call such a 5 a differential operator associated with f. Our method is 
described in full details in Section 3. Moreover, this procedure has been implemented using the 
computer algebra system Macaulay2. 

Assume that char(fc) = p > 0 and that [k : k p ] < oo. For e ^ 1 let R pe be the subring of R 
consisting of all p e - th powers of elements in R, which can also be viewed as the image of the e-th 
iteration of the Frobenius endomorphism F : R R. We set R p := R. It is shown in [Yek92, 1.4.9] 
that Vr is equal to the increasing union U e j, 0 End RP e (R). Therefore, given 6 £ Vr, there exists 
e ^ 0 such that 5 £ End rp s (R) but 5 ^ End( R ) for any e' < e. Given a non-zero polynomial 
f £ R, we have seen above that there exists 5 £ Vr that is associated with /. We say that / has 
level e if such <5 is R pe -linear, and there is no R pe -linear differential operator 5 ', with e' < e, that 
is associated with /. 

In Section 4, we study the case when / is a monomial; indeed, in Theorem 4.2 we determine the 
level of /, and we give an explicit description of the differential operator 6 associated with /. We 
also describe explicitly I e {f p _1 ), the ideal ofp e -th roots of f p -1 , where e is the level of /. The ideal 
I e (/ pe_1 ) can be defined as the unique smallest ideal J C R such that fP e ~ l £ j\P e I = (jP e | j £ J) 
(see for example [BMS08, Definition 2.2]). In Section 5 we present some families of polynomials 
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which have level one, and we give some examples. In Section 6 we focus on Elliptic Curves C C Pj| 
where F p is the finite field with p elements. We prove the following characterization: 

Theorem 1.1. Let p £ Z be a prime number and let C C P| be an elliptic curve defined by a 
cubic f(x,y,z) G F p \x,y,z\. Then 

(i) C is ordinary if and only if / has level one. 

(ii) C is supersingular if and only if / has level two. 

All computations in this article are made using the computer software Macaulay2 [GS]. 


2. Preliminaries 


The goal of this section is to review the definitions, notations and results that we use throughout 
this paper. Unless otherwise specified, k will denote a perfect field of prime characteristic p. Under 
this assumption, it is known (see [Gro67, IV, Theoreme 16.11.2]) that the ring of A;-linear differential 
operators over R = k[x\,... ,Xd\ can be expressed in the following way: 


V R \= R (D x . jt \ i = l,...,d and t > 1), where D x . t 


1 d l 
t\ dx\' 


This allows us to regard Vr as a filtered ring. Indeed, one has that 

Vr = Rr\ where := R (D Xit \ i = 1 , ..., d and 1 ^ t ^ p e — 1 ) . 

e^O 


Moreover, it is shown by A. Yekutieli (see [Yek92, 1.4.9]) that T>r = End^ P e (R), hence the previous 
filtration does not depend on the choice of coordinates. 

Now, we fix additional notation; given an a = (a±,... ,ad) € N d we shall use the following 
multi-index notation: 

x x 1 x d . 

With this notation, we set ||x Q || := max{ai,... ,a^}. By abuse of notation, we will sometimes also 
use ||a|| instead of ||x Q '|j. For any polynomial g G k[x i,..., Xd], we define 

||g|| := max ||x“||, 

x“Gsupp(g) 


where if g = Sa^- a (such that g a = 0 for all but a finite number of terms) the support of g 

is defined as 

supp(fiO := {x a G R | g Q / 0} . 

Moreover, we also define deg(g) as the total degree of g. Finally, for any ideal J C R 1 J& e 1 
will denote the ideal generated by all the p e -th powers of elements in J, or equivalently the ideal 
generated by the p e -th powers of any set of generators of J. 


2.1. The ideal of p e - th roots. Due to the central role which the ideal of p e -th roots plays 
throughout this article, we review some well-known definitions and facts (cf. [AMBL05, page 465] 
and [BMS08, Definition 2.2]). 

Definition 2.1. Given g G R, we set I e (g ) to be the smallest ideal JCU such that g G J 

Remark 2.2. Assume that k is perfect. In our assumptions, the ring R is a free f? pe -module, with 
basis given by the monomials {x“ | ||a|| ^ p e — 1}. If we write 

9 = x “’ 

0<||a||s£p e -l 

then I e (g ) is the ideal of R generated by all the elements g a [BMS08, Proposition 2.5]. 
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Remark 2.3. Notice that, if g is a homogeneous polynomial, then, for all e £ N, I e (g ) is a 
homogeneous ideal. Indeed, if we write g = ]C(K|H|^p c -i 9°* X< A then we can assume without loss 

of generality that every g(/ x" has degree equal to deg(<?). But then g a must be homogeneous of 
degree 

deg (g) - deg(x Q ) 


deg (g a ) = 


pH 


Since I e (g) is generated by the the elements g a , it is a homogeneous ideal. 

We have the following easy properties (see [AMBL05, Lemma 3.2 and Lemma 3.4] for details). 

Proposition 2.4. Given / £ R a non-zero polynomial, and given e ^ 0, the following statements 
hold. 

(i) I e (f) = I e+1 (fP). 

(ii) JeCP"- 1 ) D I e+1 {p e+1 ~ l ). 

Note that part (ii) of Proposition 2.4 produces the following decreasing chain of ideals: 


( 2 . 1 . 1 ) 


R = HP 0 - 1 ) 2 hiP- 1 ) 2 UP 2 - 1 ) 2 UP 3 - 1 ) 2 


It is shown in [AMBL05] that under our assumptions this chain stabilizes. The smallest integer 
e € N where the chain stabilizes plays a central role in this paper. We summarize the facts that we 
will need in the following theorem. See [AMBL05, Proposition 3.5, and Theorem 3.7] for details 
and proofs. 

Theorem 2.5. Let k be a perfect field of prime characteristic p, let R = k[x \,..., Xd], and let 
/ £ R \ {0}. Define 


Is -1 




= I 


e+s 


(> 


,e+ S _i 


e := inf |s ^ 1 

Then, the following assertions hold. 

(i) The chain of ideals (2.1.1) stabilizes rigidly, that is e < oo and I e _i 1 ” 1 ^ 
for any s ^ 0. 

(ii) One has 

e = min {s > 1 | 6 I s (/P*- 1 )^ 1 } , 

and e ^ deg(/). 

(iii) There exists 5 € V^ such that 5(f pB ~ 1 ) = f p£ ~ p , or equivalently such that 6(1/f) = 1 /f p . 

(iv) There is no 5' € \ with e! < e, such that 6'(1/f) = 1 /f p . 

Motivated by Theorem 2.5, we make the following definition. 

Definition 2.6. For a non-zero polynomial / £ R, we call the integer e defined in Theorem 2.5 
the level of f. Also, we will say that 6 G such that 6(f pe ~ 1 ) = f pS ~ p , or equivalently such that 
6(1/f) = l/P, is a differential operator associated with f. 

3. The algorithm 

Let A; be a computable perfect field of prime characteristic p (e.g., k is finite). Let R = 
k[x i,... ,Xd], and let / E R be a non-zero polynomial. We now describe in details the algorithm 
that computes a differential operator 6 £ Vr associated with /. 

• Step 1. Find the smallest integer e £ N such that I e (f pe ~ p ) = I e (f pB ~ 1 ). There is an 
implemented algorithm for the computation of the level of a given polynomial f £ R. Here 
follows a description: 
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Algorithm 3.1. Let k be a computable perfect field of prime characteristic p (e.g., finite), 
let R = k[x i,..., Xd], and let / £ R. These data act as the input of the procedure. Initialize 
e = 0 and flag = true. While flag has the value true , execute the following commands: 

(i) Assign to e the value e + 1, and to q the value p e . 

(ii) Compute I e (f q ~ 1 ). 

(iii) Assign to J the value I e (/ 9_1 )^. 

(iv) If f q ~ v £ J, then flag = false ; otherwise, come back to step (i). 

At the end of this method, return the pair (e, I e ( f p _1 )). Such e is exactly the e described 
in Theorem 2.5, i.e., the level of /. 


Remark 3.2. Since the level e is always at most deg(/), we can ensure that the While loop 
in Algorithm 3.1 finishes after, at most, deg(/) iterations. Notice that, a priori, there is a 
black box in this algorithm, namely the computation of / e (/ <?_1 ) at each step. However, 
the calculation of the ideal of p e -th roots is well known (cf. [KS12a, Section 6]). 

Remark 3.3. As pointed out by E. Canton in [Canl5, Definition 2.3], the so-called non- 
F-pure ideal of / introduced by O.Fujino, K.Schwede and S.Takagi in [ ST11, Definition 
14.4] turns out to be I e (Z^^ 1 ), where e is the level of / (see [FST11, Remark 16.2]). 
Therefore, Algorithm 3.1 provides a procedure to calculate the non-F-pure ideal. 

• Step 2. For e £ N as in Step 1 write / pE_1 = Ya =i Mi) where {pi, ..., p n } is the basis 

of R as an R^-module consisting of all the monomials xf 1 ■ ■ ■ x d d , with ^ p e — 1 for all 
i = 1 ,,d. 

Claim 3.4. For all i = 1,... ,n there exists <5* £ V^ such that 6i(pj) = 1 if i = j and 
5i(Pj) = 0 if i / j. 

Proof. For i £ {1,..., n} and pi = xf 1 ■ ■ ■ x a d d consider Vi := x\ 1 ai • • • x v d 1 which is 
a monomial in R because ak ^ p e — 1 for all k = 1,..., d. Notice that ViPj = (xi • • • Xd) pe ~ l 
if and only if i = j. Then set 

^i := (ii Dx k ,p e -i^ ' ,y i € ■ 

Notice that di(pii) = 1, and that if gj = x b f ■■■x b d , then Si(pj) = 0 if bk < a^ for some 
k £ {1,... ,d}. So let us assume that ak ^ b^ ^ p e — 1 for all k = l,... ,d, and that 
there is s £ {1 ,...,d} such that a s < b s . Note that by definition of ig we have that 
ViPj = x r f ■ ■ ■ x r d , with p e ^ r s ^ 2 p e — 2, so that we can write r s = p e + n for some integer 
n with 0 ^ n ^ p e — 2. 


Subclaim 3.5. The coefficient of D XstP e_i(viPj) is 


p e + n\ 

v e - V 


= 0 


mod p. 


Proof. As a consequence of a theorem proved by Lucas in [Luc78, pp. 51-52], it is enough 
to check that at least one of the digits of the base p expansion of p e — 1 is greater than the 
corresponding digit in the base p expansion of p e + n. The base p expansion of p e — 1 is 
given by 

p e - 1 = (p- 1)(1 +p-\ -hp e_1 ) = (;P ~ 1 )p° + {p~ l)p X H-1 -ijp — 1 )p e_1 , 

so that the subclaim is proved unless the first e digits of p e + n are p — 1 as well. But in 
this case, since p e + n > p e — 1 we would get 

P e + n ^ (p - 1 )p° + (p - 1 )p l H-h (p - 1 )p e ~ l +p e = 2p e - 1, 
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a contradiction since n ^ p e — 2. □ 

The Subclaim shows that D Xa = 0 for all jij with j ^ i. Therefore, using that 

Si G Vji = End/- iP e (/?,), we get 


Si{f pE X ) = ^ X) 4^ = S 4 6i ^o) = 4 ■ 

\j = 1 ) 3 =1 

(e) 

Step 3. Since 1 G T> R , for e G N as in Step 1 we have 

/P E -P G V { 4(fP e -P) = I e (F e -P) [pe] = hiP*- 1 )^ = (Cl, ■ ■ ■ ,C n )^. 


□ 


(e) 


In particular there exist a\,a n G i? such that f pE p = X^=i a 4 ■ Consider Si G V R 
as in Step 2, so that 5i(f pB 1 ) = , and set 


S:=^aiSi eV%\ 


1=1 


With this choice we have 


S(f pe x ) = S [ J24^3 = 4 ai5i ^ = 4 p , 

<3 = 1 / *J=1 


2=1 


(e) 

and using that <5 G V R we finally get 

4 (?) = 4 44 ) = 7^ 4 ^ = 


fP e ~P 1 


4. The monomial case 

Throughout this section, let k be a perfect field and let i? = /c[xi,... ,Xd\- We now analyze the 
case when f £ R is a monomial. First we show a lower bound for the level of /. 


Lemma 4.1. Let / = x^ 1 ■ ■ ■ x a d d be a monomial in R = k[x \,..., Xd), with a* > 0 for all* = 1,..., d. 
Let S G V^ R be such that S (1 //) = 1 / f p . Then, setting a := ||/|| = max {a,; | 1 ^ i ^ d}, we have 

e ^ [ lo gp(°)l + 1- 

Proof. It suffices to show that for t : = |~log p (a)] we have It(f pt p ) 2 h(f pt X ). This is because if 
the chain 

/t(/ pt - pt_1 ) 2 it(f pt ~ pt ~ 2 ) 2 ... 2 it{f pt ~ p ) 2 W^ 1 ) 
stabilizes before such step, then it would be stable at it as well, and the smallest s such that 
I s (f pS ~ p ) = l s (/ p X ) is precisely the level of / (see Theorem 2.5). Notice that t = 0 if and only if 
a% = 1 for all i. For such a monomial the Lemma is trivially true. So let us assume that t ^ 1, or 
equivalently that a* ^ 2 for at least one i G {1,..., d}. Let 

ji ■■= min{j G N | jp t ^ a*}, 


and notice that j t ^ a* for all i, and by choice of t we have that ji ^ 2 for at least one i. say j i ^ 2. 
Then 


Since (ji 

h(f pt ~ p ) 


"*-p _ p‘ai-pai . . . p t a d -pa d 


f P ~ P = X P 


= (*l 


“l-Jl 


■ ■ ■ X 


a.d~3d\p z 
d > 


■ X 


jip-pai 

1 


• - - X 


,3dP~pa d 


— 1 )p t < cii by definition of j t . we have jijf — pa* < p 4 + a* — pa* < p t . This shows that 


_ rai -31 


= Or 


x a d 3d'y Q n {.^g h an( ]; 


— 1 _ p a\—a\ 


f p - 1 = x p 


■ X 


,p t a d —a d _ f ai—1 
d 


= « 


• X 


a d 1 \pt . „P*-ai 


1 


■ X 


,p t —ad 
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which makes sense because p* ^ a* for all i. This shows that It{f pt 1 ) = (x“ x 1 • • • x a d d x ), and 
because j\ ^ 2 we have 


It(f pt ~ p ) 


(a? 1- * 


rr a d jd 
X d 


)D(x 


ai — 
1 


2 


“d- 1 

X d 


) 2 (* 


ai — 1 
1 


“d-i 

x d 


) = Jt(/ P *- 1 ). 


□ 


Theorem 4.2. Let / = x^ 1 • • • x^ d be a monomial in fc[xi,..., Xd\, with > 0 for all i = 1,..., d. 
Let a := \\f\\ = max{aj | 1 ^ i ^ d}. Then / has level e := |~log p (a)~| + 1, and I e (f pe ~ 1 ) = 
(x? 1 ” 1 • • • x“ d_1 ). Furthermore, 


*=-n 


P^-pCLi 


‘ Ar;,p e -1 ' x i l 


ai—1 


2—1 


is a differential operator associated with /. 


<E X> 


( e ) 

R 


Proof. Set e := |~log p (a)] + 1. During the proof of Lemma 4.1 we already proved that I e (/ pe 1 ) = 
(x“ 1_1 • • ■ x a d d ~ 1 )] keeping this fact in mind, it is enough to check that 5(f pe ~ 1 ) = fP e ~P. Indeed, 
we have 


SifP*- 1 ) = 5 (xf 01-01 


x a 


p e a d -a d \ _ ( x ai~l 


■ X 


a d -l 


5 f x : 


„P ~a\ 


JP ~a d \ _ 


n 

2—1 


X P ai-pat 


D^pe^xf ^ 


_ p a\ —pa i 


= X 


1 


• X‘ 


p c a d -pa d _ fp e -p 


and therefore the proof is completed. 


□ 


Regarding the level e obtained in Theorem 4.2, one might ask whether, given any non-zero f £ R, 
its level would always be bounded above by |~log p (||/||)] + 1. Unfortunately, this is not the case, 
as the following example illustrates. 

Example 4.3. Consider / := xy 3 + x 3 £ F 2 [x,y]. In this case, one can check with Macaulay2 
that the level of / is 4, while |"log 2 (||/||)l +1 = 3. In fact, the level is even strictly greater than 
[log 2 (deg(/))] + 1 = 3. 

The monomial in Theorem 4.2 is assumed to be of the form x“ x • • -x^ d . Using a suitable linear 
change of coordinates, we immediately get the following Corollary, which includes the general 
monomial case. 


Corollary 4.4. Let n ^ d, let / = if 1 ■ ■ ■ be a product of powers of linear forms which are 
linearly independent over k, and let a := max{aj | 1 ^ i ^ n}. Then / has level e = |~log p (a)~| + 1, 
the ideal of p e -th roots is I e {f pe ~ l ) = and 

n 

»n ■ du*-i ■ <?-*) 6 v r 

2—1 

is a differential operator associated with /. Here, if li = ^ ij x ji then the differential operator 

Rt, ,p e —i is defined as ff j i Ay ^xj ,p e — i • 


5. Families of level one 


Polynomials of level one, that is polynomials / such that I\{f p ~ 
For instance, let f,g £ R and let 5 £ T>r be associated with /. 
5' := 5 ( g p^ 1 ■ _) we get 


1 ) = R , are somehow special. 
Assume that e = 1, then for 


6' 





p 
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The authors do not know whether, for any choice of /, g £ R, f 0, there always exists S' £ Dr such 
that S'(g/f) = (g/f) p - In fact, when 6 G for e ^ 2, the best we can get is S' (g/f) = g pE // p , 
with S' := 5 (g pe_1 • _). On the other hand, for any / G R we have Rf = Dr • j and therefore, for 
any g G R, there exists S' G Dr such that S' {1/ f) = g p / f p . In fact it is enough to set S' := g p ■ 5. 

We will now exhibit some families of polynomials that have level one, together with some exam¬ 
ples. However, before doing so, we want to single out the following elementary statement, because 
we will be using it repeatedly throughout this section. It may be regarded as a straightforward 
sufficient condition which ensures that a polynomial has level one. In this section, unless otherwise 
stated, k will denote a perfect field and R = k[x i,..., Xd] will be a polynomial ring over k. 

Lemma 5.1. Let f € R be a non-zero polynomial, write 





and assume that fp is a unit for some 0 ^ ||/3|| ^ p e — 1 and some e ^ 1. Then, / has level one. 

Proof. By definition, we have that I e {f pe i ) = R ; on the other hand, we know that I\{f p ~ l ) D 
I e (f p 1 )- In this way, combining these two facts it follows that I\{f p ~ l ) = R , and therefore / has 
level one. □ 


We can give an easy but useful characterization of homogeneous polynomials of level one. 

d 

Lemma 5.2. Let f £ R be a homogeneous non-zero polynomial. Let {pj}j = \ ■= {x Q | ||a|| ^ p— 1} 
be the monomial basis of R as a A p -module. Then / has level one if and only if pj G supp(/ p_1 ) 
for some j = 1,... ,p d . 


Proof. Write 




E ft x ° 

0^||a||<p— 1 


Note that, since / is homogeneous, Ii(/ P “ 1 ) is a homogeneous ideal by Remark 2.3. If / has 
level one, then Ii(f p ~ l ) = (f a \ 0 ^ ||a|| ^ p — 1) = R, therefore there exists at least one 
coefficient fp that is outside of the irrelevant maximal ideal m = (xi,..., xf)- Write fp = X + r, 
with A G k and r G m. Then, we can write / p ~ l = A p x^ + h for some h G R. Also, since 
{x“ | 0 ^ ||a|| ^ p — 1} is a basis of R as R p -module, there is no cancellation between A p x^ and h. 
Thus, pj = x' 3 G supp(/ p ^ 1 ). Conversely, if pj = x /3 G supp(/ p_1 ) for some 0 ^ ||/3|| ^ p — 1, then 
we can write again / p_1 = A p x 3 + h for some A G k and some h G R. Also, we can assume that 
x i:i ^ supp (h). Then the coefficient f p of x' 3 in the expansion of / p_1 must be A p + r p = (A + r) p for 
some r G m, and thus A + r G /i(/ p_1 ). Since the latter is homogeneous (here we are using again 
Remark 2.3), we have in particular that A G ii(/ p_1 ), which implies /i(/ p ^ 1 ) = R, and therefore 
/ has level one. □ 


Proposition 5.3. Let f £ R be a non-zero polynomial whose support contains a squarefree term 
involving a variable that does not appear in any other term of the support of /. Then / has level 
one. 


Proof. Without loss of generality we can assume that x\---x n £ supp(/), and that x\ does not 
appear in any other term of supp(/). Write / = Aoaq • • • x n + where mi = x ai and 

ati £ N d are of the form (0, ,..., a^) for * = 1,..., s. Then, 

P - 1 = E ( ■ P ~ 1 ■ ) A o° Al 1 • • • Ah*! ■ ■ ■ XnPmf ■■■ml*. 

i 0 +ii+~+is=p-l ’ ’ ’ ls ' 
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Notice that in order for a term in the support of / p—1 to be divisible by it is necessary that 
*o = p — 1, in which case i\ = ... = i s = 0. Hence x\~ l ■ ■ ■ Xn -1 is in the support of / p_1 , appearing 
with coefficient Aq 1 ^ 0. Then we can write 

f p ~ 1 = Ar 1 *r 1 ---zr 1 + E 

0 < ||a|| ^ p - 1 
ai ^ p - 1 

and the Proposition now follows from Lemma 5.1, since Aq 1 is a unit. □ 

Example 5.4. Let / = x 2 + y 2 + xyz £ ¥ p [x,y, z]. Since £ appears in the square free term xyz 
and nowhere else in the support of / we have that / has level one by Proposition 5.3. In fact, 
5 := D X 'p-iD y ^p_iD Z) p-i £ is a differential operator associated with /. 

Proposition 5.5. Let f £ R be a non-zero polynomial of degree n such that every element of its 
support is a squarefree monomial. Then / has level one. 

Proof. Without loss of generality we can assume that x\ - ■ -x n £ supp(/). We want to show that 
we can apply Lemma 5.1. Let / = Aooq • • • x n + ^* =1 ^i m £ where the monomials nq are squarefree 
of degrees di := deg(mi) ^ n, that we can assume different from x\ ■ ■ ■ x n . Then 

f- 1 = E (■ P ~ 1 ■ ) A o°-M 1 ''' (*! ‘ ‘ ‘ xnrml 1 ■ ■ ■ . 

... • . VO) ■ ■ ■ ) *s/ 

*0+*iH- Hs=p— 1 

Note that the choice io = p — 1, i\ = ... = i s = 0 gives the monomial Aq 1 {x\ ■ ■ ■ x n ) p_1 , and we 
want to show that this choice of indices is the only one that gives such a monomial. By way of 
contradiction, assume that (aq • • • x n ) lo rnf • • • m l s s = (aq • • • x n ) p_1 , then necessarily each rrn divides 
aq • • • x n , because they are squarefree. Since we are assuming that none of the monomials mi is 
equal to aq • • • x n , we must have that deg(mj) < n. But then 

deg ^(aq • • ■ x n ) l °mf ■ ■ ■ m 1 /^ = ni 0 + d\i\ H-b d s i s < n(i 0 + i\ 4-f if) = n(p - 1), 

which is a contradiction because (aq ■ ■ ■ x n y°m l f ■ ■ ■ rn l s s = (aq ■■■x n ) p ~ 1 , and the degree of the 
latter is n(p — 1). Therefore if we write 

/£x“, 

0^||«||=£p—1 

then the coefficient of 1 • • • 1 is precisely Aq - 1 , which is a unit. Using Lemma 5.1, the 

Proposition now follows. □ 

Example 5.6. Let R = F p [Xy]i^i be a polynomial ring in n 2 variables and let / = det(Wj)- 
Because of Proposition 5.5 / has level one, since its support consists only of squarefree monomials. 

Proposition 5.7. Let f £ R = k[x\,...,Xd] be a homogeneous quadric. Then / has level one 
unless / is the square of a linear form, in which case / has level two. 

Proof. If / is a power of a linear form, then / has level two by Corollary 4.4. Otherwise, if p 2 
there exists a linear change of variables that diagonalizes / (cf. [Ser73, Chapter IV, Proposition 
5]). That is, we can assume that, after a linear change of coordinates, / = x 2 + • • • + ax 2 , where 
2 ^ n ^ d and a is either 1 or an element of k which is not a square. Notice that x^ l x^ 1 appears 

with coefficient A := (p_i) £ k \ {0} if n ^ 3, and with coefficient q(p -1 )/ 2 ( p _i) £ k \ {0} if n = 2. 

2 2 
Therefore / has level one by Lemma 5.2. Finally, if p = 2 and / is not a power of a linear form, 

then we can assume that X 1 X 2 appears with non-zero coefficient in / p_1 = /, and we conclude using 
again Lemma 5.2. □ 
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Proposition 5.8. Let / = x\ + • ■ • + x^ £ R be a diagonal hypersurface of degree t ^ 1. 
t ^ min{d,p} and p = 1 (mod t), then / has level one. 


Proof. Our assumptions on t,d and p allow us to expand f p 


-l 


f p ~ L = 


(P !)• p— l p— l 


p-i 

t 


i 


+ 


E 


in the following manner: 
(P 1)" ti\ tij. 


*lH-h«d=p-l 




If 


Since ,^ p _ 1 )\t £ k \ {0}, the above equality shows that x p ■ ■ ■ x p 1 appears in supp(/ p 
\\~r~ )■) 

non-zero coefficient, hence / has level one by Lemma 5.2. 


— In 


with 

□ 


The assumptions p = 1 (mod t) and t ^ min{d, p} in Proposition 5.8 cannot be removed in 
general, as the following examples illustrate. 

Example 5 . 9 . Let R := ¥^,[x,y, z] and / := x 3 + y 3 + z 3 . One can check using Macaulay2 that / 
has level two. Notice that, in this case 3 = deg(/) ^ min{3, 5} and 5 = 2 (mod 3). 

On the other hand, consider now R : = Fy[x, y\ and / := x 3 + y 3 . One can check using Macaulay2 
that / has level two. In this case 3 = deg(/) > 2 = min{2, 7} and p = 7 = 1 (mod 3). 

The diagonal hypersurface considered in Proposition 5.8 is of the form x\ + • • • + using a 
suitable linear change of coordinates, we immediately get the following Corollary, which includes 
as a particular case Proposition 5.8. 

Corollary 5 . 10 . Let n ^ d, let / = £\ + • • • + £f n be a diagonal hypersurface of degree t ^ 1 made 
up by linear forms £±,... ,£ n which are linearly independent over the field k. If t ^ min{n, p} and 
p = 1 (mod t), then / has level one. 

Before going on, we want to review the following notion (see [ILL + 07, page 243]): 


Definition 5.11. A polynomial f £ R is said to be regular provided 

w): = 

where Tj(/) denotes the Tjurina ideal attached to /. 


In characteristic zero, a polynomial is regular if and only if its Bernstein-Sato polynomial is 
bf(s) = s + 1 [ L + 07, Theorem 23.12], In this case, Rf is generated by 1// as a D-module. 

Proposition 5.12. Let fc be a perfect field of characteristic 2, and let / £ k[x i,..., xj be regular. 
Then, / has level one. 


Proof. Since / is regular, there are ro, ri,..., £ R such that 

df 


1 = r 0 j +n- -h 

OX I 


+ r d 


dx d 


In this way, setting 


it follows that 5(f) 


d 

5 -= r o + ^2 

3 =1 

1 and therefore / has level one. 


d 

dxj 


□ 


Remark 5.13. A very easy way to produce polynomials which are simultaneously regular and 
of level one in arbitrary prime characteristic works as follows. Let fc be a perfect field of prime 
characteristic p , let R = k[x i,..., x^], and assume that f £ R is a non-zero polynomial of the form 
/ = Ax* + g , for some 1 ^ i ^ d, some A £ k — {0}, and some g £ R such that 

vp- = 0 (i.e., g £ k[x i,...,Xj_i,x),x i+ i,x i+2 ,... ,x d }). 
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Then, / is regular and of level one; indeed, the fact that / is of level one follows directly from 
Proposition 5.3. 


6. Elliptic Curves 

Let p £ Z be a prime and let C C be an elliptic curve defined by an homogeneous cubic 
f(x,y,z ) £ F p [x,y,z\. We want to review here the following notion (see [Hus87, 13.3.1]). 

Definition 6 . 1 . C is said to be ordinary if the monomial ( xyz ) p_1 appears in the expansion of 
fP~ l with non-zero coefficient. Otherwise, C is said to be super singular. 

The general form of a cubic defining an elliptic curve is the following 

/ = y z + a\xyz + a$yz — x — a, 2 X z — a^xz — a%z , 

where a \,..., clq £ F p . When p / 2,3, the expression above can be further simplified to 

/ = y 2 z - x 3 + axz 2 + bz 3 , 

for a, b £ F p (see [Hus87, 3.3.6] for details). We are now interested in computing the level of elliptic 
curves C. We are mainly interested in upper bounds, since it is easy to see from Lemma 5.2 that any 
ordinary elliptic curve has level one, and that any supersingular elliptic curve has level at least two. 
First, we explore the low characteristic cases, where the list of possibilities (up to isomorphism) is 
very limited. 

Proposition 6 . 2 . Let C C P| be a supersingular elliptic curve defined by a cubic / £ F p [x,y,z\. 
If p = 2 or p = 3, then / has level two. 

Proof. Set D := D xp 2 _^D yp 2 _iD zp 2 _^. By [IIus87, 13.3.2 and 13.3.3] there are only the following 
two cases, up to isomorphism: 


p 

Elliptic curve 

Differential operator 

2 

x 3 + y 2 z + yz 2 

y 2 Dx 3 z + z 2 Dx 3 y + x 2 Dxyz 2 

3 

S 2 2 

x — xz — y z 

(x e z 3 — x 3 y 6 )Dx 4: z 5 + 

+(x 9 + x 3 z 6 + y 6 z 3 )Dxy 8 + y 3 z 6 Dx 4 y 5 


The table above is exhibiting a differential operator of level two for each polynomial, therefore the 
level is at most two in all such cases. We have already noticed that C is ordinary if and only if / 
has level one. This shows that when p = 2 or p = 3, C is supersingular if and only if / has level 
two. □ 

Recall that R = F p [x, y, z] is a free /?, p2 -module with basis given by {x r y s z t \ 0 ^ r, s, t ^ p 2 — 1}. 
For a polynomial g £ F p [x, y, z] consider 

9= (c(r,s,t)) p2 x r y s z\ 

and recall that by Proposition 2.4, with this notation, 12 (g) is the ideal generated by the elements 
c(r, s, t ), for 0 ^ r, s, t ^ p 2 — 1. 

Remark 6.3. If / £ F p [x, y, z] is a cubic and g = / p2_1 , then for any 0 ^ r, s, t ^ p 2 — 1 one has 


deg (f p2 x ) 


3 (p 2 ~ 1) 

p2 


p2 


deg(c(r, s, t)) ^ 


= 2 . 
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In particular, l 2 {f p2 *) = (c(r, s,t) | 0 ^ r, s, t ^ p 2 — l) is generated in degree at most two. 

For the rest of the section, we will denote ii(/ p_1 ) and / 2 (/ p2 “ 1 ) simply by I\ and 1 2 . 

6.1. Preliminary computations. The purpose of this part is to single out some technical facts 
which will be used for proving the main result of this section; namely, Theorem 6.9. 

Lemma 6.4. Let / = y 2 z — x 3 + axz 2 + bz 3 £ F p [x, y, z], where p ^ 2, 3. Then 

c(0,p 2 - 2 ,p 2 - 1) = y. 

Proof. A monomial in the expansion of f p2 ~ l will have the form {y 2 z) h (—x 3 ) 1 (axz 2 y (bz 3 ) k , where 
h + i + j + k = p 2 — 1. Looking at the coefficient of y p ~ 2 z p in such expansion, by degree 
considerations we only have three possibilities: 

{ ; X P 2 . yP 2 ~ 2 z p2 ~ 1 

y p 2 . yP 2 ~ 2 Z P 2 - ] 

Z P 2 . yP 2 -2 Z P 2 -1 

Since p 2 — 2 is not even, there is no choice of h that realizes the first and the third cases. So we 
are left with the second, which is achieved only by the choice h = p 2 — 1, i = j = k = 0. This 
shows that the coefficient of y p ~ 2 z p ~ l in the expansion of f p ” 1 is precisely y p , and the Lemma 
follows. □ 


Before going on, we need to review the following classical result, due to Legendre, because it 
will play some role later in this section (see Proof of Lemma 6.7). We refer to [AZ04, page 8] for a 
proof. 

Theorem 6.5 (Legendre). Let n ^ 0 be a non-negative integer, let p be a prime number, and let 
(jp{n ) be the sum of the base p digits of n. Then, 

/ n _ n ~ °>( n ) 


p t | m}. 


p- 1 

where, given any non-negative integer m ^ 0, 

v p (m) := max{f ^ 0 : 


Lemma 6.6. Let p / 2 be a prime. Then, 

(p 2 - 1)! 


A := 


H-ill 


/o 


(mod p). 


Proof. On one hand, p 2 — 
since p / 2 it follows that 


1 = (p — 1)(1 + p) is the base p expansion of p 2 



(¥■ 


(1 +p) 


1; on the other hand, 


is the base p expansion of (p 2 — l)/2. Keeping in mind these two facts it follows, using Legendre’s 
Theorem, that 


Vp (A) = v p {fp 2 


1)!) - 2 v p 



p 2 -l-2(p-l) 

P-1 



= P -1 - (p -1) = 0 , 


hence p does not divide A and therefore we can ensure that A / 0 (mod p). 


□ 
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Lemma 6.7. Let p / 2,3 be a prime, and let / = y 2 z — x 3 G F p [x, y, z]. Then, I\ = I 2 = (x,y). 
In particular, / has level two. 


Proof. Consider the expansion 


p-1 


r - 1 = E - 


2—0 


(P r * r 3(p 2 -l-*) 

i!(p2-l-i)! J ' 


For i = ( p 2 — l)/2 we obtain the monomial 

A/-^(P 3 -!)/ 2 ^-!)/ 2 = A x p2 • (x^- 3 )/ 2 /- 1 ^- 1 )/ 2 ) , 

where A = ^1, 2 ^ 0 (indeed, this follows by Lemma 6.6). Because of the term z 1 in the 

r) 

expansion above, the choice i = (p 2 — l)/2 is clearly the only one that gives the monomial 
x (p 2 -3)/2yp 2 -i 2; (p 2 -i)/2 q £ £j ie p as j s Q f pi as an /7,P 2 - mo dule. Therefore c f 3 ,p 2 — 1. = 

A 1 / p2 x = Ax, and thus x G I 2 . In addition, by Lemma 6.4 we always have y & 1 2 . Therefore 
(x,y) C I 2 . On the other hand, consider the expansion 

f P-i = 


p-i 

^ j!(p 2 - 1 - j)\ 


^ P y 2j z j x 3 ^- l ~ j) . 


We claim that either 2 j ^ p or 3(p— 1— j) ^ p. In fact, suppose j < p/2, or equivalently j ^ (p—1)/2, 
because j is an integer. Then 3(p — 1 — j) ^ 3(p — l)/2 ^ p since p ^ 5 by assumption. This 
shows that all the coefficients c(r, s, t) p in the expansion of / p_1 = Y^o^r s t^p-i c ( r > s > t) p x r y s z t are 
contained in (x,y)^, and thus /1 = (c(r, s,t) | 0 ^ r,s,t ^ p — 1) C (x,y). Therefore the Lemma 
follows from the chain of inclusions (x,y) C I 2 C Ii C (x,y). □ 

Lemma 6.8. Let p / 2,3 be a prime and let C be a supersingular elliptic curve defined by 
/(x, y, z) = y 2 z — x 3 + axz 2 + 6z 3 . If either a = 0 or b = 0, then / has level two. 

Proof. Notice that, since C is supersingular, we have F p [x, y, z] = R D I\ D J 2 , and we want to show 
that I\ = I 2 . Also, by Lemma 6.4 we always have y G I 2 . If a = b = 0, then Lemma 6.7 ensures 
that I\ = I 2 = (x,y). Now assume a / 0 and 6 = 0. We claim that c f 1 ,p 2 — 1, = pz, 

2-1 


•j 2 z) h (—x 3 ) 1 (r r "" 2 ^ 
L z~— the possibilities are 


for some g / 0. A monomial in the expansion of f p 1 will have the form (y 2 z) h ‘{—x 3 ) l {axz 2 y for 

9 p 2 - 1 2_1 »> 2 —'i 

h + i + j = p — 1. Looking at the terms involving x 2 y p 

„2 

X‘ 

y 2h x 3i+j z 2j+h = 


The second case is not possible, since 2p 2 — 1 is not even and hence here is no h that gives it. This 
forces h = (p 2 — l)/2, and hence the first case is also not possible, since for such an h the exponent 

of z must be at least (p 2 — l)/2. Canceling y p2 ~ 1 z S ~^~ we are left with 


P 2 " 1 „2. 

2 Q 

1 P-3 

X 2 yP 

2) 2 

P 2 - 1 „2. 

2 0 

1 P -3 

X 2 yP 

2) 2 

P 2 " 1 p 2_ 

2 Q 

1 P -3 

X 2 yP 

L Z 2 


x 3 ^‘ = x^zP 2 - 1 


r.2_ 1 . 


This implies that j = p 9 1 and i = 0. Therefore the coefficient ofx 2 / 1 z 2 in the expansion 
of / p2_1 is 

=a^z p2 = (pzf, 
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where fi = I a 2 




7 ^ 0. This shows the claim. With similar considerations, one can see that 

P 2 -3 2 P 2 ~ 

- 7 T- ,P ~ 1 ,-A- I = X. 


This gives (x, y, z) C. I 2 C. I\ C R 1 which forces I\ = I 2 because I\ is a proper homogeneous ideal 
of R = ¥ p [x,y,z\, hence I\ C ( x,y,z ). Finally, if a = 0 and 6 7 ^ 0 the arguments are completely 
analogous to the case a / 0, 6 = 0. Here we get 

„2 Q J2 


P - 3 2 , P ~ 1 

,P - 1, 


= x 


and c ( 0 ,p 2 — l,p 2 — 2 ) = Az 


2 ’2 

for some A 7^ 0. Therefore ( x,y,z ) C / 2 C T C R. which once again forces I\ = I 2 . 
6.2. Main result. Next statement is the main result of this section. 


□ 


Theorem 6.9. Let p £ Z be a prime number and let C C P|^ be an elliptic curve defined by a 
cubic f(x,y,z ) € F p \x,y,z\. Then 

(i) C is ordinary if and only if / has level one. 

(ii) C is supersingular if and only if / has level two. 

Proof. Part (i) follows from Lemma 5.2, which also shows that if / has level at least two, then C 
is supersingular. So it is left to show that if C is supersingular, then / has level precisely equal to 
two. By Proposition 6.2 we only have to consider the cases where p 7 ^ 2, 3, thus we can assume 
that / is of the form 

f(x, y, z) = y 2 z — x 3 + axz 2 + bz 3 . 

for a,b G ¥ p . Furthermore, by Lemma 6.8 we can assume that ab 7 ^ 0. First we claim that 


p 2 - 3 

2 


,P 2 - 1, 


v-1) - 


'j = x + A z, where 




(p 2 ~ 1 )! 


(-1 )*a^“ 3i 6 2i+1 g ¥ p 


In fact a general monomial in the expansion of f p2 1 will have the form (y 2 z) h (—x 3 ) 1 (axz 2 y (bz 3 ) k , 


T^o )!*! *¥-3i !(2i + l)! 


and by looking at terms that involve x 2 y- 
possibilities: 


2h x 3i+j z h+2j+3k _ < 


p 2 -l 


Z 2 


by degree considerations we have three 


2 0 

P -3 

p 2 -i ^ 

1 2 

2/P 2 

2 Q 

P -3 

„2 1 ^ 

to 

2/P 2 

2 Q 

P —3 

p 2 -i 

2 

2/P 2 2 


The second case cannot be realized since 2 p 2 — 1 is not even, hence we necessarily have h = 
This leaves two cases: 


p 2 -i 


x 3i+j z 2j+3k _ 


X 2 


• X 2 


The first one can only happen when h = i = p » and j = k = 0, giving the monomial 


/ 9 ' V - 1 . Q. P -! , , . P - 1 „2 p--J t) 2_1 P~-J- „2_1 

( 2 / 2 ) 2 (-X ) 2 = (-1) 2 X r . X 2 y P L z 2 = X P ■ X 2 yf L z 2 
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because p 2 = 1 (mod 4) . For the second case i, j and k will have to satisfy 

J = £!_3 - 3z 

k = 2i + 1 

As both j and k must be non-negative, we have 

p 2 — 3 p 2 — 7 


0 ^ i € 


6 


6 ’ 


because p 2 = 1 (mod 6). Furthermore, the coefficient of (y 2 z) h (—x 3 y(axz 2 y(bz 3 ) k under these 
conditions is precisely 


(p 2 - 1)! 


(^ 1 ) ! * ! (^-3i) ! (2t + l)! 


(_i y a p 2 3 3 i 6 2 *+i, 


proving that 


p 2 -3 2 n 2 -l xP ' 


->p - V 


^ = x p2 + Xz p2 = {x + X 1 ^ 2 z) p2 


for A as above. Since we are working over F p , we finally have A 1/,p2 = A, showing the claim. With 
similar arguments, one can see that 


p 2 — 1 2 o p 2 +l\ i 

- 3, ) =ax + fiz, 


c (^r-iP 2 — i) = ~ ax + TZ -> 

p 2 — 1 

c (0,p 2 — l,p 2 — 2) = Ox + 6“ 2 ~ 


where 0 G F. 






(p 2 - 1 )! 


and 


r = 


U (£^ 3 )!i!(£^i-3z)!(2i + l)! 


« (p3 - 1)! 

(^i)m(^-3i)!(2z + 3)! 


(-l)*a^~ 3 V m € F p , 


(-1 ) i a 2L ^' 3i 6 2i+3 G Fp. 


Since /2 = (c(r, s, t) | 0 ^ r, s, t ^ p 2 — 1), in particular we have 

P 2 - 1 

(x + Az, ax + /uz, — ax + rz, 0x + 6 2 z) C I 2 . 

Claim. The following matrix has full rank: 

' a fi 

1 A 

—a t 

p 2 - 

0 b~ 
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In fact if det 


a n 
1 A 


7 ^ 0 then we are done, otherwise we have p = aX. Note that 

P 2 r i 

C v 2 - 1)! 


P = 


E 


p^—3 

i =0 


!i!(^i- 3 z)!( 2 z + l)! 


(-l) i a E r 1 - 3i 6 2< + 1 


E 


(p 2 - 1 )! 


jo 2 —3 

i=0 1-2“ 


!z! -3*)!(2* + l)! 

_(p 2 ~ 1)!_ 

i=o l ^)h ! (^i-3z)!(2f + l)! 




, f (P 2 -!)! , 

1 / p 2 -3\i /v-iy /v+ 2 \, 


P -7 
6 

E 


(-1 )*a £ ^- 3 i 6 2i+1 


3(p 2 - 1)! b ^+ 

^ / P 2 1 A i ( P 2 "7 \ | ( p 2 +2 A , 


where the last equality comes from the fact that p 2 = 1 (mod 12 ) and by rearranging the binomial 
coefficients. Also, 


« A = E 


(p 2 - 1 )! 


i =0 




Using that p = aA, we get 


2^-3ij!(2z + l)! 
3(p 2 - 1)! 




P 2 ~M | / P 2 ~7 \ | / P 2 +2 \ , 


b 3 


E 


1=1 L 2 


1 1 


p 2 — 1 p 2 —1 


- 3i 


(P 2 - 1)! 




3i ) \{2i + 1)! 


-(-l) i a E -2 1 - 3i b 2i+1 


E 


—3 (p 2 — 1)! 


i (^)!(i-l)!(A 1 “3i)!(2i + l)! 


(-1 ya^-H 2 ^ 1 


3 V -f— —r— 

( p 2 - 1A i .• i 


(p 2 - 1 )! 


(-lya^-^b 2 ^ 3 , 


* (^i)h!(2^-3z)!(2f + 3)! 
where the last equality comes from reindexing the sum. Since 3 is invertible in F p we get that 


0 = 


E 


(p 2 - 1)! 


(—l) 2 a P 2 3*522+3 


*=o l 2 




3z) !(2z + 3)! 


(p 2 — 1) ! P 2 +2 


p 2 -i y fp 2 -7y / p a +2 a , 


b 3 
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p _7 

because p 2 = 1 (mod 12), hence (—1) — e~ = —1. But the latter is precisely r, and this argument 
shows that if oA = p, then t = 0. Therefore either aA / p or 


—a t 


—a 0 

p 2 -1 

= det 

p 2 - 1 

0 b *-3- 


6 b V- 



Hence the matrix has rank two, and the Claim follows. 

p 2 -i 

But this shows that there are linear combinations of x + A z, ax + pz, —ax + tz and Ox + b 2 z 
that produce x and z, that is ( x,z ) C I 2 . By Lemma 6.4 we always have that y € 12- Therefore 


(■ x,y,z ) C I 2 C Ii C R, 


implying that I\ = I 2 and hence that / has level two. 


□ 


7. A Macaulay2 session 

The purpose of this section is to explain, through a Macaulay2 session, how the algorithm 
introduced in Section 3 works in specific examples. 

We begin clearing the previous input and loading our scripts. 

clearAll; 

load "differentialOperator.m 2 "; 

We fix the polynomial ring that we will use throughout the following examples. 

p= 2 ; 

F=ZZ/p; 

R=F [x,y,z,w]; 

The first example illustrates a particular case of Theorem 4.2. 
i 6 : f=x~3*y~5*z"7*w"4; 


i7 : L=differentialOperatorLevel(f); 


i 8 : L 


2 4 6 3 

08 = (4, ideal(x y z w )) 

This means that, in this case, / has level 4 and that I^f 2 ^^ 1 ) = (x 2 y 4: z 6 w 3 ). Now, we produce a 
differential operator 6 of level 4 such that 5(1/f) = 1 /f 2 . 

i7 : DifferentialOperator(f) 

o7 = | xl0y6z2w8 d_0~15d_l~15d_2~15d_3~15x2y4z6w3 | 

As the reader will note, the output is a row matrix; it means that 5 turns out to be 
x 10 y 6 z 2 w 8 ■ D x ^_ 1 D y)2 i_ l D Z)2 i_ l D w ^-i ■ x 2 y 4 z 6 w 3 . 

Our next aim is to illustrate a particular case of Corollary 4.4. 

11 8 : f=x~3*(x+y)~5*(x+y+z)~7*(x+y+z+w)”4; 

119 : L=differentialOperatorLevel(f); 
iilO : first L 


00 IO = 4 
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Now, a particular case of Proposition 5.3. 

1113 : f=x~2+y~2+z~3+x*y*z*w; 

1114 : L=differentialOperatorLevel(f); 

1115 : L 

ool5 = (1, ideal 1) 

This means that / has level one. Now, we produce the corresponding differential operator. 

i7 : DifferentialOperator(f) 

o7 = | 1 d_0d_ld_2d_3 | 

It means that the differential operator produced in this case is D Xt \D y ^D z ^D Wt i. 

The following computation may be regarded as a particular case of Proposition 5.5. 

16 : f=x*w-y*z; 

17 : DifferentialOperator(f) 
o7 = | 1 d_0d_ld_2d_3yz I 

It means that, in this case, the differential operator produced is D x \D y \D z ]D w \ ■ yz. 
Next, a homogeneous quadric (cf. Proposition 5.7). 

16 : f=x~2+y~2+x*y+z~2+w~2; 

17 : DifferentialOperator(f) 

o7 = | 1 d_0d_ld_2d_3zw | 

We finish with a homogeneous cubic. 

1112 : f=x~3+y~3+z~3+v~3; 

1113 : L=differentialOperatorLevel(f); 

1114 : L 

ool4 = (2, ideal (w, z, y, x)) 
ool4 : Sequence 

1115 : DifferentialOperator(f) 

ool5 = | w2 d_0~3d_l"3d_2~3d_3~3x3z3w | 

I z2 d_0"3d_l~3d_2"3d_3"3x3zw3 | 

I y2 d_0~3d_l~3d_2"3d_3~3yz3w3 | 

I x2 d_0~3d_l~3d_2"3d_3~3xy3z3 I 

This means that our differential operator in this case turns out to be 

w 2 ■ D x , 3 D y}3 D z>3 D w>3 ■ x 3 z 3 w + z 2 • D x , 3 D y 3 D z 3 D V! 3 ■ x 3 zw 3 
+ y 2 ■ D Xt3 Dy )3 D Zt3 D W)3 ■ yz 3 w 3 + x 2 ■ D x , 3 D y - i D z ^D vx3 ■ xy 3 z 3 . 
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8. The code of the algorithm 

The aim of this final section is to show our implementation in Macaulay2 of the algorithm 
described in Section 3 of this manuscript; the whole code can be found in [ 5DSV14]. Throughout 
this section, R = F p [xi,..., will be the polynomial ring with d variables with coefficients in the 
field Fp. 

First of all, we write down the code of a procedure which, given an ideal I of R, return as output 
/[ pe l, i.e., the ideal generated by all the p e th powers of elements in /. The method below is based 
on code written by M. Katzman and included, among other places, in [ \S12b]. 

frobeniusPower(Ideal, ZZ) := (I,e) ->( 

R:=ring I; 
p:=char R; 
local u; 
local answer; 

G:=first entries gens I; 
if (#G==0) then 
{ 

answer=ideal(0_R); 

> 

else 

{ 

N:=p~e; 

answer=ideal(apply(G, u->u"N)); 

>; 

answer 

); 

Now, we exhibit the code of a function which, given ideals A, B of R, produces as output the ideal 
I e {A) + B. For our purposes in this manuscript, B = (0) and A is a principal ideal. Once again, 
this is based on code written by Katzman and included in [KS12b]. 

ethRoot(Ideal,Ideal, ZZ) := (A,B,e) ->( 

R:=ring(A); 
pp:=char(R); 

F:=coefficientRing(R); 
n:=rank source vars(R); 
vv:=first entries vars(R); 

Rl:=F[vv, Y_l. .Y_n, Monomial0rder=>Product0rder-[n,n}-,MonomialSize=>16] ; 
q:=pp~e; 

JO:=apply(1..n, i->Y_i-substitute(vv#(i-l)~q,Rl)); 

S:=toList apply(l..n, i->Y_i=>substitute(vv#(i-l),R1)); 

G:=first entries compress( (gens substitute(A,Rl))%gens(ideal(JO)) ); 

L:=ideal 0_R1; 
apply(G, t-> 

{ 

L=L+ideal((coefficients(t,Variables=>vv))#1); 

»; 

Ll:=L+substitute(B,Rl); 

L2:=mingens Ll; 

L3:=first entries L2; 

L4:=apply(L3, t->substitute(t,S)); 
use(R); 
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substitute(ideal L4,R) 

); 

Next, we provide the code of our implementation of Algorithm 3.1. Namely, given f £ R, the 
procedure below gives as output the pair (e, I e ( f p ~ 1 )), where e is the level of /, and I e ( f p _1 ) is 
the ideal where the chain (2.1.1) stabilizes. As the reader can easily point out, this method is just 
turning Theorem 2.5 into an algorithm. 

differentialOperatorLevel(RingElement):=(f) ->( 

R:=ring(f); 
p:=char(R); 
local J; 
local I; 
local e; 
e=0; 

flag:=true; 
local q; 
local N; 
while (flag) do 
{ 

e=e+l; 
q=p~e; 

N=q-1; 

I=ethRoot(ideal(f~N),ideal(0_R),e); 

J=frobeniusPower(I,e); 

N=q-p; 

if ((f"N)7„ J==0) then flag=false; 

>; 

(e ,1) 

); 

Now, let x = Xi (for some 1 ^ i ^ d), n ^ 0, and / € R. The below procedure returns as output 

1 d n f 
n! dx n 

It is worth noting that, in some intermediate step of this method, we have to lift our data to 
characteristic zero in order to avoid problems with the calculation of 1/n!. 

DiffOperator(RingElement,ZZ,RingElement):=(el,numb,funct)-> 

( 

R:=ring(el); 

vv:=first entries vars(R); 

S:=QQ [vv]; 

functl:=substitute(funct,S); 
el=substitute(el,S); 
for i to numb-1 do 
functl=diff(el,functl); 
functl=l/(numb!)*functl; 
funct=substitute(functl,R); 
el=substitute(el,R); 
use R; 

return funct; 

); 
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Next, we provide a method which, given a monomial x Q = X T ' ‘ ‘ x °d 0 ^ a* ^ p - 1 for all i, 

returns as output the differential operator 5 a £ T>^ where, for any other monomial = x\ J 1 • • • x^ 
with 0 ^ bi ^ p — 1 for all i, 5 a acts in the following way: 

j„(*/>) = I 1 ’ if “ = A 

V / I 0, otherwise. 

As the reader can easily point out, the below method is just turning Claim 3.4 into an algorithm. 
DeltaOperator(RingElement,ZZ):=(el,pe)-> 

( 

R:=ring(el); 
indet:=vars(R); 

for i to numColumns(indet)-1 do 

( 

deg:=degree(indet_(0,i),el); 
listdiff_i=pe-l-deg; 

); 

return listdiff; 

); 

The next two methods are quite technical; however, both are necessary in order to avoid problems 
during the execution of our main procedure, which we are almost ready to describe. 
checkcondition=method(); 

- checkcondition finds a random monomial in startingpol//pol with all the 

- variables in degree <p~e 

checkcondition(ZZ,RingElement,Ideal,ZZ):=(pe,startingpol,J,indexvar) -> 

( 

pol:=(first entries gens J)_indexvar; 
genJ=first entries gens J; 

R:=ring(pol); 
var:=vars R; 
nvars:=numColumns var; 
for i to numgens J-l do 

sett_i=set first entries monomials(startingpol//genJ_i); 
supportpol= first entries monomials (startingpol//pol); 
contat=0; 

for i to #supportpol-l do 

( 

flag=true; 

for k to nvars-1 do 

if degree(var_(0,k),supportpol_i)>pe-l then flag=false; 
if flag then 
( 

rightsupport_contat=supportpol_i; 
contat=contat+l; 

); 

); 

correctmon=rightsupport_(random contat); 
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); 


return correctmon; 


DifferentialAction=method(); 

- DifferentialAction computes differential(element) 

DifferentialAction(RingElement,RingElement):=(differential,element)-> 

( 

R:=ring(element); 

T:=ring(differential); 
var:=vars T; 
nvars:=numColumns var; 
total=0; 

diffmonomials=first entries monomials (differential); 
for i to #diffmonomials-1 do 
( 

moltiplication=substitute(element,T); 
for j to floor(nvars/2)-l do 

moltiplication=moltiplication* 

var_(0,floor(nvars/2)+j)~(degree(var_(0.floor(nvars/2)+j),diffmonomials_i)); 
differentiation=coefficient(diffmonomials_i.differential)*moltiplication; 
for j to floor(nvars/2)-l do 

differentiation=DiffOperator(var_(0,nvars-1-j), 
degree(var_(0,floor(nvars/2)-l-j),diffmonomials_i).differentiation); 
total=total+differentiation; 

); 

return total; 

); 

We conclude showing our implementation of the algorithm described in Section 3, which is the 
main result of this paper. 

DifferentialOperator(RingElement):=(f)-> 

( 

R:=ring(f); 
p:=char(R); 

F:=coefficientRing(R); 
local e; 
local J; 

(e,J)=differentialOperatorLevel(f); 

J=frobeniusPower(J,e); 

variable:=vars(R); 

variablel:=first entries variable; 

genJ:=first entries gens J; 

nvars:=numColumns vars R; 

powerf:=f"(p"e-l); 

T:=F[d_0..d_(nvars-1).variablel]; —creating ring of differentials 
varDelta:=first entries vars T; 
powerfT=substitute(powerf,T); 
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use R; 

condition:=true; 
while condition do 
( 

for i to numgens J-l do 

( 

listsupport=checkcondition(p"e,powerf,J,i); 
listdiff=DeltaOperator(listsupport,p~e); 
deltal=l; 
delta2=l; 

for j to nvars-1 do 

( 

use T; 

deltal=deltal*varDelta_j ~(p"e-l); 
expon2=listdiff_j; 

delta2=delta2*varDelta_(nvars+j)"expon2; 

); 

delta_i=deltal*delta2; 

newgen_i=DifferentialAction(delta_i,powerfT); 

); 

newgenmatr=matrix(newgen_0); 
for i from 1 to numgens J-l do 

newgenmatr=newgenmatrInewgen_i; 
newgenmatrR=substitute(newgenmatr,R); 
if J==ideal(newgenmatrR) then 
( 

matalpha=f"(p"e-p)//newgenmatrR; 
matalpha=substitute(matalpha,T); 
for i to numgens J-l do 

if i==0 then matrixdiff=matrix-[-[matalpha_(i ,0) ,delta_i}]- 
else matrixdiff =matrixdiff | |matrix{-[matalpha_(i,0) ,delta_i]-]-; 
totale=0; 

for i to numgens J-l do 

totale=totale+matalpha_(i,0)*DifferentialAction(delta_i,powerf); 
powerfp=substitute(f“(p"e-p),T); 
return matrixdiff; 
condition=false; 

); 

); 

use R; 


Acknowledgements 

This project started during the summer school PRAGMATIC 2014. The authors would like 
to thank Aldo Conca, Srikanth Iyengar, and Anurag Singh, for giving very interesting lectures 
and for sharing open problems. They also wish to thank Alfio Ragusa and all the organizers of 
PRAGMATIC 2014 for giving them the opportunity to attend the school. The authors also thank 
Josep Alvarez Montaner and, once again, Anurag Singh and Srikanth Iyengar, for many helpful 
suggestions concerning the material of this manuscript. 


AN ALGORITHM FOR DIFFERENTIAL OPERATORS IN POSITIVE CHARACTERISTIC 


23 


References 

[AMBL05] J. Alvarez Montaner, M. Blickle, and G. Lyubeznik. Generators of D-modules in positive characteristic. 
Math. Res. Lett, 12(4):459-473, 2005. 1, 2, 3 

[AZ04] M. Aigner and G. M. Ziegler. Proofs from The Book. Springer-Verlag, Berlin, third edition, 2004. Including 
illustrations by Karl H. Hofmann. 11 

[BDSV14] A.F. Boix, A. De Stefani, and D. Vanzo. differentialoperator.m2: a Macaulay2 package for 
computing a differential operator which acts as the Frobenius homomorphism. Available at 
http://atlas.mat.ub.edu/personals/aboix/research.html, 2014. 18 

I. N. Bernstein. Analytic continuation of generalized functions with respect to a parameter. Funkcional. 
Anal, i Prilozen., 6(4):26-40, 1972. 1 

M. Blickle, M. Mustafa, and K. E. Smith. Discreteness and rationality of E-thresholds. Michigan Math. 

J. , 57:43-61, 2008. Special volume in honor of Melvin Hochster. 1, 2 

R. Bpgvad. Some results on D-modules on Borel varieties in characteristic p > 0. J. Algebra, 173(3):638- 
667, 1995. 1 

E. Canton. A note on injectivity of Frobenius on local cohomology of hypersurfaces. Available at 
http://arxiv.org/pdf/1502.03184vl, 2015. 4 

O. Fujino, K. Schwede, and S. Takagi. Supplements to non-lc ideal sheaves. In Higher dimensional algebraic 
geometry, RIMS Kokyuroku Bessatsu, B24, pages 1-46. Res. Inst. Math. Sci. (RIMS), Kyoto, 2011. 4 
A. Grothendieck. Elements de geometrie algebrique. IV. Etude locale des schemas et des morphismes de 
schemas IV. Inst. Hautes Etudes Sci. Publ. Math., (32), 1967. 2 

D. R. Grayson and M. E. Stillman. Macaulay2, a software system for research in algebraic geometry. 
Available at http://www.math.uiuc.edu/Macaulay2/. 2 

D. Husemoller. Elliptic curves, volume 111 of Grad. Texts in Math. Springer-Verlag, New York, 1987. 
With an appendix by Ruth Lawrence. 10 

S. B. Iyengar, G. J. Leuschke, A. Leykin, C. Miller, E. Miller, A. K. Singh, and U. Walther. Twenty-four 
hours of local cohomology, volume 87 of Grad. Stud. Math. American Mathematical Society, Providence, 

RI, 2007. 1, 9 

M. Katzman and K. Schwede. An algorithm for computing compatibly Frobenius split subvarieties. J. 
Symbolic Comput., 47(8):996-1008, 2012. 4 

M. Katzman and K. Schwede. FSplitting, a Macaulay2 package implementing an algorithm for computing 
compatibly Frobenius split subvarieties. Available at http://katzman.staff.shef.ac.uk/FSplitting/, 2012. 

18 

E. Lucas. Sur les congruences des nombres euleriens et les coefficients differentiels des functions 
trigonometriques suivant un module premier. Bull. Soc. Math. France, 6:49-54, 1878. 4 
J.-P. Serre. A course in arithmetic, volume 7 of Grad. Texts in Math. Springer-Verlag, New York- 
Heidelberg, 1973. 8 

A. Yekutieli. An explicit construction of the Grothendieck residue complex. Asterisque, (208):127, 1992. 
With an appendix by Pramathanath Sastry. 1, 2 

Department of Economics and Business, Universitat Pompeu Fabra, Jaume I Building, Ramon Trias 
Fargas 25-27, 08005 Barcelona, Spain. 

E-mail address', alberto.fernandezb@upf.edu 
URL : http://atlas.mat.ub.edu/personals/aboix/ 

Department of Mathematics, University of Virginia, 141 Cabell Drive, Kerchof Hall Charlottesville, 
VA 22903, USA 

E-mail address', ad9fa@viginia.edu 

URL', http: //people . virginia.edu/~ad9fa/ 

Departimento di Matematica e Informatica, Universita di Firenze, Viale Morgagni, 67/a - 50134 
Firenze, Italy 

E-mail address : davide.vanzo@unifi.it 


[Ber72] 

[BMS08] 

[B0g95] 

[Canl5] 

[FST11] 

[Gro67] 

[GS] 

[Hus87] 

[ILL+07] 

[KS12a] 

[KS12b] 

[Luc78] 

[Ser73] 

[Yek92] 


